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Denote by H{t) ^ {Hi{t), . . . , HN{t)) a function in f £ R+ with values in (0,1)'^. Let 
{B"'-^\t)} = {B"'-*\t),te R+ } be an (A^, d)-multifractional Brownian sheet (mffls) with Hurst 
functional H{t). Under some regularity conditions on the function H{t), we prove the existence, 
joint continuity and the Holder regularity of the local times of {B^^*^ (t)}. We also determine 
the Hausdorff dimensions of the level sets of {_B^'*'(f)}. Our results extend the corresponding 
results for fractional Brownian sheets and multifractional Brownian motion to multifractional 
Brownian sheets. 

Keywords: Hausdorff dimension; level sets; local times; multifractional Brownian sheets; 
one-sided sectorial local non-determinism 

1. Introduction 

A one-dimensional fractional Brownian motion (fBm) ^" = {^°'{t),t E R} with Hurst 
index a € (0, 1) is a real-valued, centered Gaussian process with covariance function 
given by 

E[r(s)rw] = |[kp" + itp"-it-sp"] vs,tGM. (1.1) 

It was introduced, as a moving average Gaussian process, by Mandelbrot and Van Ness 
[31]. 

Fractional Brownian motion has interesting properties such as self-similarity of or- 
der a € (0, 1), stationary increments and long-range dependence (when a > 1/2), which 
make it a good candidate for modeling different phenomena in, for example, finance 
and telecommunication. However, this model may be restrictive, due to the fact that all 
of its regularity and fractal properties are governed by the single Hurst parameter a. 
To model phenomena whose regularity evolves in time, for example, Internet traffic or 
image processing, Levy-Vehel and Peltier [28] and Benassi, Jaffard and Roux [7] indepen- 
dently introduced multifractional Brownian motion (mfBm) in terms of moving average 
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representation and harmonizable representation, respectively. Multifractional Brownian 
motion is governed by a Hurst functional a{t) with certain regularity in place of the 
constant Hurst parameter a in fBm. 

Several authors have studied sample path and statistical properties of mfBm. For 
example, Benassi, Jaffard and Roux [7] considered the sample path Holder regularity of 
mfBm and determined the Hausdorff dimension of its graph. Ayache, Cohen and Levy- 
Vehcl [2] and Hcrbin [22] studied the covariancc structure of mfBm from its harmonisable 
representation. Recently, Boufoussi, Dozzi and Gucrbaz [13, 14] studied the existence, 
joint continuity and Holder regularity of the local time of mfBm and established Chung's 
law of the iterated logarithm for mfBm. The main tool they applied to derive their results 
is the property of one-sided local non-dcterminism of mfBm. 

There are multiparameter extensions of fBm, among which, two typical examples are 
multiparameter Levy fBm ^" = {£,°'{t),t £ M^} (whose covariancc function is given by 
(1.1) with I • I being the Euclidean norm in R^) and fractional Brownian sheets, where the 
former is isotropic, while the latter are anisotropic in general. Since they were introduced 
by Kamont [24] (see also Ayache, Leger and Pontier [4]), fractional Brownian sheets (fBs) 
have been studied extensively as a representative of anisotropic Gaussian random fields 
in recent years. See, for example, Dunker [18], Mason and Shi [32], 0ksendal and Zhang 
[33], Xiao and Zhang [45], Ayache and Xiao [6], Ayache, Wu and Xiao [5], Wu and Xiao 
[39] and the references therein for further information. Still, the regularity of fBs does 
not evolve in the TV-dimensional 'time' parameter t S M.{^. 

To model anisotropic Gaussian random fields whose regularity evolves in time, such as 
images, Ayache and Leger [3] and Hcrbin [22] introduced so-called multifractional Brow- 
nian sheets (mfBs) in terms of their moving average representations and harmonisable 
representations, where the constant Hurst vector of fBs is substituted by Hurst function- 
al. Furthermore, they showed that mfBs have a continuous modification and determined 
the pointwise and local Holder exponent of mfBs. They also proved that mfBs are locally 
self-similar. We refer to Ayache and Leger [3] and Hcrbin [22] for the definitions of the 
corresponding concepts and results. 

In studying anisotropic random fields, Xiao [44] suggested that it is more convenient 
to use the following metric p^^, on : 

N 

p^{s,t)=J2\si-t,\''' Vs,teM^, (1.2) 
e=i 

where K=iKi,...,KN) € (0,1)^ is a fixed vector. Denote by H{t) = {Hi{t), . . . , HN{t)) 
a function in t e M.^ with values in (0, 1)^. We say that H{t) satisfies Condition A if 
there exist a positive number a G (0, 1) and a vector {Ki, . . . , K^) S (0, 1)^ such that 

A.l. for every i € {1, . . . , N}, a < Hi{t) < Ke for aU t e R^; 

A. 2. Hi{t) {i = l,. . . ,N) satisfies a /9^,-Lipschitz condition on every compact set, that 
is, for every compact subset / C R;'^, there exist positive constants = (/) and 
S such that 

\Hi{t) - Hi{s)\ <c,p^{s,t) Vs,i e / with \s-t\< S. 
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Now, we are ready to define multifractional Brownian sheets via their moving average 
representations. 

Definition 1.1. Let H{t) = {Hi{t), . . . , HN{t)) be a function in t G M.^ with val- 
ues in (0, 1)^ satisfying Condition A. A real-valued multifractional Brownian sheet 
{B^^*\t)} = {B^^*\t),t g R^} with functional Hurst index H(t) is defined as the fol- 
lowing moving average Wiener integral: 

r ^ 

Bo^'\t)= r\[{te-ue)f^''>~'/^~{~u,)'^'^'^~'/^]W{du) VteR^, (1.3) 

where s+ = max{s, 0}, and where W — {W{s), s G M^} is a standard real-valued Brown- 
ian sheet. 

Remark 1.2. Our Definition 1.1 generalizes the definition in Ayache and Leger [3], 
where they define Hi{-) as a function of ti G M+ (£ = 1, . . . , N). Herbin [22] defines multi- 
fractional Brownian sheets by using the following moving average representation (cf. [22] 
Definition 2, page 1261): 

r ^ 

B'o"'-'\t)^ rr[|if-w,|^^W-i/2-|i.,|^"(*)-i/2]M/(dii) VieM^. (1.4) 
Jr" 

Based on Dobric and Ojeda [16] (see also Stoev and Taqqu [35]), we know that the multi- 
fractional Brownian sheet defined by (1.3) and the multifractional Brownian sheet defined 
by (1.4) have different correlation structures in general, even if iV= 1. Our definition is 
more convenient to use than that of Herbin's when we derive the one-sided sectorial local 
non-determinism for multifractional Brownian sheets in Section 2. The form (1.3) is pre- 
ferred in some applications in the one-dimensional case because it is easier to separate 
the future from the past. 

Define the Gaussian random field {B^(*)(t)} = {B^(*)(i) :t e E^} with values in R'* by 

B"^'\t) = {B"^'\t),...,B"'^'\t)) yteR^, (1.5) 

where {Sf (t)}, . ■ . , {sf (t)} are d independent copies of {B(f (t)} . Then (t), 
tG M.^} is called an (A^, d) -multifractional Brownian sheet with fmictional Hurst index 
Hit). 

Note that if A^ = 1, then {B^'^*^^{t)} is a multifractional Brownian motion in M"^ with 
Hurst index Hi{t) £ (0, 1); if A^ > 1 and Hi{t) = Hi, . . . , HN{t) = Hn , then {B^(*)(<)} is 
an (A^, (i)-fractional Brownian sheet with Hurst index H = {Hi, . . ., Hjsi), which will be 
denoted {B" {t),t eR^}. 

In this paper, we will study the existence and regularity of the local times of multifrac- 
tional Brownian sheets. Our main technical tool is the property of one-sided sectorial local 



868 



M. Meerschaert, D. Wu and Y. Xiao 



non- determinism. This property is more general than the sectorial local non-determinism 
which was first introduced by Khoshnevisan and Xiao [27] for the Brownian sheet and 
then extended by Wu and Xiao [39] to fractional Brownian sheets. See Section 2 for more 
information. 

Our results show that multifractional Brownian sheets are similar to fractional Brow- 
nian sheets in many ways. They admit jointly continuous local times when the indices 
Hi{t) stay in the range for which the fractional Brownian sheet has jointly continuous 
local times. We also establish a Hausdorff dimension result and, essentially, the dimension 
of the level set is the same as for the constant parameter case, except that we take the 
supremum of the constant parameter formula. We also show that the supremum can be 
taken locally, to establish that the fractal dimension of the random field varies in space. 
Hence, multifractional Brownian sheets are useful in applications such as composite ma- 
terials, or porous media flow, when the material properties vary in space. They may also 
find useful applications in image processing. 

The rest of this paper is organized as follows. In Section 2, we prove some basic results 
on mfBs that will be useful to our arguments. In Section 3, we provide a sufficient condi- 
tion for the existence of L^-local times of the (iV, d) mfBs and prove that the condition 
also implies the joint continuity of the local times. We prove the Holder regularity of the 
local times in Section 4. Finally, we derive the local Hausdorff dimensions of the level 
sets of {B^^*\ty\ in Section 5. Our results extend the results of Ayache and Xiao [6] 
and Ayache, Wu and Xiao [5] for fractional Brownian sheets and Boufoussi, Dozzi and 
Guerbaz [13, 14] for multifractional Brownian motion to multifractional Brownian sheets. 

We end the Introduction with some notation. Throughout this paper, the underlying 
parameter space is K.^ or = [0, oo)^. A parameter t e is written a,s t = {ti, ... ,tN), 
or as (c) , if ti = • • • = t^r = c. For any s, t € such that Sj < tj [j = 1, . . . , N) (wc denote 
this by s ^ t) , we define the closed interval (or rectangle) [s, t] = Il^i [^3 ' ^j] ■ We use A to 
represent the collection of closed intervals [s,i] with € [£,A]^ for some fixed positive 
numbers e and A. For any integer m > 1, we always write A„i for the Lebesgue measure 
on R™ and use (•,•) and | • | to denote the ordinary scalar product and the Euclidean 
norm in M™, respectively. 

Throughout this paper, an unspecified positive and finite constant will be denoted by 
c, which may not be the same in each occurrence. More specific constants in Section i 
are numbered as c. ^ , c. ^ , . . . . 

2. Preliminaries 

In this section, we will provide some lemmas that are useful for proving our main results. 
Lemma 2.1 is an extension of Boufoussi et al. [13], Lemma 3.1 from fractional Brownian 
motion to fractional Brownian sheets. 

Lemma 2.1. Let < £ < A and 0<a<7<l he fixed constants. Let {Z^{t), {t,K) G 
M.^ X [a, 7]^} be a real-valued Gaussian random field defined by equation (1.3) with 
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H{t) = K. There then exists a constant c, ^ = c(a, 7, e, A, TV) > such that 

E[Z^it)~zi^'it)f<c,JP-P'\' (2.1) 
for allte [e,A]^ and all f3,(3' e [a, 7]^. 

Proof. For any /3,/3' G [a, 7]^, wc define — (3 and — {I3[,... . . . ,/3Ar) for 

j = 1, . . . , TV. Clearly, = /3'. Since 

N 

nZ^^{t) - Z^\t)r < NY,E[Zr\t) ZS\t)]\ (2.2) 

it suffices for us to prove that for j S {1, . . . , N} fixed, 

E[Zo"'" (t) - ZS' {t)f < c{a, 7, e, A, N)\f3, -P'/. (2.3) 
By using the moving average representation of fBs (equation (1.3)), we have that 

E[zr\t)-zs'{t)r 

i-i 

.l3'e-l/2.2 



AT 

X n [{U-w)l^-'/'-i-u,)l^-'/'rdu, (2.4) 

X / [it, - .,)r^/^ - i-u,)r'^' - ((., - .,)5-^/^ - hu,)r'^Y 

Jr 
■ =h X 

where iij = (ui, . . . , Uj+i, . . . , ujv)- 
A change of variables shows that 

i-i N 

^i=^n*f'- n (2-5) 

which is bounded for all t G [e, A]^ and /3, ^9' e [a, 7]^. 

Next, we estimate l2- Without loss of generality, we may assume f3j > (i'y We rewrite 
I2 as the following summation: 



'\it,-u,f'-'/^-{t,~u^^'.-'/^Y^, 



00 
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(2.6) 



:=IIi+II: 



2- 



By the mean value theorem, we have that for some P'j < rjji < f3j [rjji may depend on 



III = r [{h - Mt, - u.j)f\P, - d.., 

" (2.7) 

<c,,|/3,-/3;|2 

for all t E [£,A]^ and /3,/3' G [a, 7]^. In the above, the last inequality follows from a 
change of variables. 

Similarly, we have that for some f3j < rjj2 < Pj, 

112 - f [{t, - u,r^-^l^ ln(i, - u,) - (-u,)"^-^-^/' ln(-%)]^|/3, - An, 

(2.8) 

for all t e [e,A]^ and e [0,7]^. Equation (2.3) is proved by combining (2.4)-(2.8). 
This proves Lemma 2.1. □ 

Combining Lemma 2.1 and Lemma 8 of Ayache and Xiao [6], we have the following 
result. 

Lemma 2.2. Let {Sq (i)} he a multifractional Brownian sheet in M. There exist pos- 
itive constants 6 >0, ^ and ^ such that for all s, < G [e, A]^ with \s — t\ < S, for any 



u 



S rifci l^i ^ti,Si\/ tt] , we have 



N N 

e=i 1=1 
Proof. By the elementary inequalities 

3(a2 + fe^ + c^) > (a + + cf > - Ah^ - 4c^ 

we have that 

3(E[So^("^ it) - Bo""^"' {s)f + ¥.[B^^'^ {t) - B^^""^ {t)f + E[B(f (s) - B^^""^ {s)f) 
>E[B"^'\t) ^ B"^'\s)f 



(2.10) 

-4E[Bf"'(.)-i?f^'(.)]l 



> iE[B^("^ (t) - {s)f - 4E[B^(*' {t) - {t)f 
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The first term on the right-hand side of (2.10) is the variance of the increment of a 
fractional Brownian sheet with Hurst index H(u). By Lemma 8 in Ayache and Xiao [6], 
there exist positive constants c^ g and , depending only on a, 7 = maxjiCi, . . . , Kn}, e 
and N ^ such that 

N N 

c.,« E - '^\'"''-^^ ^ E[i?„^(")(t) - B^'^-\s)f < c,, \t, - (2.11) 
1=1 1=1 

Meanwhile, by Condition A and Lemma 2.1, there exists (5 > small such that for all 
t € [e, A]^ with |s — t| < 5, we have |u — s| < 5, which implies 

E[So^("' {s) - {s)f < c,,, \Hiu) - Hisr 

(2.12) 

AT N 



S, 



^,J2mu)~He{s)\^<c,,J2\t,- 



1=1 1=1 

and, similarly, 

N 

e=i 

Combining (2.10) with (2.11), (2.12) and (2.13), and noting that sup^^ Hi{ti) < Ki (still 
by Definition 1.1), we can see that there exists i5 > 0, which depends only on a,7,e, A 
and N, such that for |s — t| < S, (2.9) holds. This completes the proof of Lemma 2.2. □ 

Remark 2.3. For any fixed t g [e, A]^ and any fixed unit vector ry S R^, let us consider 
the increment 

Y^{t)^B^^'+'^\t + r^)-B^^'\t). 

In light of Lemma 2.2, it would be interesting to compute the bounds on the correlation 
function 

..,(.) = JM±^1ML (2.14) 

VE[y„(t + s)2]E[r,(t)2] 

and its time- and direction-varying spectral density (or Wigner-Ville distribution) 

/ c-'<^^«Vt,,(s)cls (2.15) 



since those are the properties of the random field that are usually used to model natural 
phenomena. 



As we mentioned in the Introduction, the main technical tool of this paper is the prop- 
erty of sectorial local non-determinism. Hence, it will be helpful to say a few words about 
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the history of various forms of local non-determinism and briefly recall some of the defi- 
nitions. The concept of local non- determinism, which provides a powerful tool for dealing 
with the complex dependence structure of Gaussian processes, was first introduced by 
Berman [11] to unify and extend his methods for studying local times of real- valued 
Gaussian processes, and then extended by Pitt [34] to Gaussian random fields. The no- 
tion of strong local non- determinism was later introduced by Cuzick and DuPreez [15] 
for Gaussian processes (i.e., A'^ = 1) and can be extended in various ways to Gaussian 
random fields (see Xiao [43]). The simplest is the following. Let X = {X{t),tG R^} be 
a real-valued Gaussian random field with < E[X(t)^] < oo for all t E I, where / e ^ is 
an interval. Let </) be a given function such that 0(0) = and (j){r) > for r > 0. Then 
X is said to be strongly locally (jj -non- deterministic (SL0ND) on / if there exist positive 
constants n, and ro such that for alH g / and all < r < min{|t|,ro}, 

Var(X(t)|X(s) ■.seI,r<\s^t\<ro)> c,,,„0(r). (2.16) 

Here, Var(X(t)|X(s), s £ S) denotes the conditional variance of X{t), given X{s), s G S. 
Pitt [34] proved that fractional Brownian motion ^" is SL(/)ND with <j>{r) =r^". When 
iV = 1 and Va.i{X {t)\X{s) s G I ,r < t — s < ro) > -^^(f){r), X is said to satisfy one-sided 
strong local (j)-non-determinism. The properties of strong local non-determinism have 
played important roles in studying the regularity of local times, small ball probabilities, 
uniform Hausdorff dimension results and other sample path properties of Gaussian pro- 
cesses and Gaussian random fields. We refer to Xiao [43, 44] for more information on 
applications of SL^ND, as well as other forms of local non-determinism. 

It is well known that the Brownian sheet does not satisfy the property of local non- 
determinism in the sense of Berman and Pitt (hence, it does not satisfy (2.16) either). 
Recently, Khoshnevisan and Xiao [27] introduced the concept of sectorial local non- 
determinism and proved that the Brownian sheet has the property of sectorial local non- 
determinism. Wu and Xiao [39] extended the result of Khoshnevisan and Xiao [27] and 
proved that the fractional Brownian sheet {B^ {t),t S where H = {Hi, . . . ,H]y) G 

(0,1), also has the property of sectorial local non-determinism. Namely, for any fixed 
positive number e G (0, 1), there exists a positive constant c^ u, depending only on e,H 
and N, such that for all positive integers n>2, and alH-'^, . . . , i" G [e, oo)^ , we have 

N 

Var(<(r)|<(ti),...,<(f-i))>c,,,,^ min (2.17) 

J = l 

where i° = 0. 

The concept of sectorial local non-dcterminism opens the door to unifying the previ- 
ously different treatments for the Brownian sheet and fractional Brownian motion, and 
can also be used to prove new results about fractional Brownian motion. See, for exam- 
ple, Khoshnevisan. Wu and Xiao [26], Wu and Xiao [39] and Ayache, Wu and Xiao [5] for 
studies on sample path properties of the Brownian sheet and fractional Brownian sheets 
by using the sectorial local non-determinism property. 
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It is not known whether the multifractional Brownian sheet {B^^*\t)} is sectorially 
locally non-deterministic. To study the local times of (A^, d)-mfBs {B^^*\t)} , we will 
work with multifractional Liouville sheets (mfLs) at first and prove that they can be well 
approximated by Gaussian processes which have the property of one-sided strong local 
non-determinism (see (2.22) and (2.23) below). For a function H{t) = {Hi{t), . . . ,iJjv(t)) 

Hit) 

satisfying Condition A, the real- valued, centered Gaussian random field {Xq ^ '(t)} = 
{X"^'\t),te R^} defined by 

f ^ 

X"^'\t)^ r\ite-se)"'^''>-'/'W{ds), teR^, (2.18) 
J[o.t] fJl 

is called a multifractional Liouville sheet with functional Hurst index H{t). One param- 
eter mfLs were first introduced by Lim and Muniandy [30] as an extension of fBm; see 
Lim [29] for more properties on one parameter mfLs. 



It follows from (1.3) that for every t E 



r ^ 

B^^'\t)=X^^'\t)+ l[gdte,sc)W{ds), (2.19) 

J(~oo,t]\[0,t] i^-^ 



where 

9i{U, s,) = {{t, - s,Uf'^'^-'/' ((-.s,)+)^^«-V2 



and the two fields on the right-hand side of (2.19) are independent. We will show that 
in studying the regularity properties of the local times of {_Bq (i)}, the Liouville sheet 
{Xq {t)} plays a crucial role and the second field in (2.19) can be neglected. More 
precisely, wc will make use of the following property. For all integers n > 2, , . . . , e R'^ 
and ui , . . . , u„ g R, we have 

Var u,B^^''\tn^ > Var (j^ UjX^^''Ht=)^ . (2.20) 

Next, we use an argument in Ayache and Xiao [6] to provide a useful decomposition 
for {X(^^*^(i)}. For every t G [e, A]^, we decompose the rectangle [0,t] into the following 
disjoint union of subrectangles: 

TV 

[0, t] = [0, U y R,{t) U A(£, t), (2.21) 

£=1 

where Ri{t) = {r G [0, A]^ -.0 < ri < e ii i ^ £,e < re < t^} and A(£, t) can be written as a 
union of 2^ — iV — 1 subrectangles of [0,t]. Denote the integrand in (2.18) by g{t,r). It 
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follows from (2.21) that for every t€ [e,oo)^, 

N 

g{t,r)W{dT) 



X^^'\t)= [ g{t,r)W{AT)+y^ [ g{t,r)W{dr) 

J[0,e]" 1=1-' Rdt) 



A(e,t) 



N 



■.= X{e,t) + Y.Ye{t) + Z{e,t). (2.22) 
1=1 

Since the processes X(e,t), Yi{t) {l<i< N) and Z{e,t) are defined by the stochastic 
integrals with respect to W over disjoint sets, they are independent Gaussian random 
fields. 

The following lemma shows that every process Yi (t) has the property of strong local 
non-determinism along the ^th direction. It will be essential to our proofs. 

Lemma 2.4. Let I £ A and let £ & {1,2, . . .,N} be fixed. For all integers n>2 and 
, . . . e I such that 

if < if < • • • < , 

we have 

Var(r,(r)|y,(i^') : < j < n - 1) > c, ,, \t^ - t^-' \^"^^'"\ (2.23) 
where t^ = and > is a constant depending only on e and I. 

Proof. The proof is in the same spirit as the proof of Lemma 2.1 of Ayache, Wu and Xiao 
[5]. Working in the Hilbert space setting, the conditional variance in (2.23) is the square 
of the L^(P)-distance of Ye{t^) from the subspace generated by Yi{t^){0 < j < n — 1). 
Hence, it is sufficient to show that there exists a constant c^ -^^ such that 

E||y,(r) - a,r.(t^) j > c,.,ji^ - tr'^^'^^ (2.24) 

for all Oj € M (j = 1, . . . , n — 1). However, by splitting Rilf^) into two disjoint parts and 
using the independence, we derive that 



/i?f(t")\flf(t"-i) 



git\r)Widr)] 



E(Y,{t'')-J2a,Ye{t^)^ >e|^^ 

> r... r ... rn(i^-r..)^^^(*")-Mr (2.25) 

Jo Jt"^^ Jo 



— ''2,12 '-e 



This proves (2.24) and hence Lemma 2.4. □ 
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Combining Lemma 2.2 and Lemma 2.4 with the proofs of Lemma 2.1 and Lemma 8.1 
in Berman [11] (see also Boufoussi, Dozzi and Guerbaz [13], Theorem 3.3), we have the 
following. 

Proposition 2.5. For every integer n > 2, there exist positive constants d and d (both 
of these may depend on n) such that for every £ = I, . . . , N , 

(n \ n 

J2uj[Mt3) - Yeit^-^)] > CnJ2^^j^iiT[Yi{t^) - Mt^-')] (2.26) 
i=i / i=i 

for all (mi, . . . , Un) G M" and all points t^ , . . . ,t" I satisfying t^ < ■ ■ ■ <t" with t" — tl< 
6. 

H(t) 

The following lemma relates the multifractional Brownian sheet {Bq {t)} to the 
independent Gaussian random fields (£ = l,...,iV), which is a direct extension of 
Ayache, Wu and Xiao [5], Lemma 2.2. 

Lemma 2.6. Let I ^ A. For all integers n > 2, t^, . . . , £ / and mi, . . . , m„ G M, we have 

Cn \ N / n \ 

Y,uX^*'\t') > Y^V^AY^u.Yeit^) . (2.27) 



Consequently, for all k G {1,...,A'^} and positive numbers pi,...,pk > 1 satisfying 
= 1, we have 



1 " 

— <n 

11/2 J-1 



[detCov(i?o^('')(ti), . . .,<(*")(t"))]V2 - l\ [detCoviYiit^),...,Y,itr^))]m^P^) ' 

(2.28) 

where det Cov(.Zi, . . . , Z„) denotes the determinant of the covariance matrix of the Gaus- 
sian random vector (Zi, . . . , Z„). 

Remark 2. 7. By (2.27) and Lemma 2.4, we derive that {B^^*^ (t)} satisfies the following 
property: for all integers n>2 and t^ , . . . ,t^ G I such that P -< for every j <n — 1, 

N 

Y^r{B^r )\B^{t'), . . . , B^{t-')) > c,,,, J] min {t] - t^)2^.(*"). (2.29) 

J = l 

This property is weaker than (2.17) and will be referred to as one-sided sectorial local 
non- determinism. 



We will also make use of the following technical lemmas, among which Lemma 2.8 is 
from Xiao and Zhang [45], Lemma 2.9 is proved in Ayache and Xiao [6] and Lemma 2.10 
and Lemma 2.11 are from Ayache, Wu and Xiao [5]. 
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Lemma 2.8. Let <h <1 be a constant. Then for any S > 2h, M > and P >0, there 
exists a positive and finite constant depending only on S,s, (3 and M, such that for 
allQ<a<M, 

-1 fi 

dr [a+\s-r\^'']-^ds<c,,,ia-^^-'/^^ + 1). (2.30) 



(2.31) 



Lemma 2.9. Let a, 13 and rj be positive constants. For A> and B > 0, let 

dt 



J:^J{A,B)= [ 
Jo 



(A + i")'3(B + t)'7' 



There then exist finite constants and depending only on aj3 and rj, such that 

the following hold for all reals A,B>0 satisfying A^^°' < c^ i^B : 



(i) if a(3 > 1, then 



(ii) if a(3 = 1, then 



- ^2,16 ^/3_a-i^^ ; {2.32) 



J<c,,,,-^log(l+i3A-i/a). (2.33) 
(iii) if < ceP < 1 and a(3 + rj^l, then 

J<c...[-B^^^+l)■ (2.34) 

Lemma 2.10. Let (p^,...,p„) C (0,1)^. For any qe [0,J2f^j^p-^), let t € {1, . . . , N} 

be the integer such that 

j:-<q<j:-, (2.35) 



e=i 1=1 P' 



with the convention that X]£=i ^ ■= 0- There then exists a positive constant A,- < 1, de- 
pending only on (p^ , . . . , p„), such that for every A G (0, A^), we can find t real numbers 
Pe>^ <t) satisfying the following properties: 

M<i V£=l,...,r (2.36) 

and 



(l-A)E— ^^^'^ + --E-- (2.37) 



1=1 7^1 
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Furthermore, if we let := X]I=i ^7 Q > 0, then for any positive number p € (0, 

there exists an Iq € ■ ■ ■ ,t} such that 

— + 2p,„P<l. (2.38) 

Lemma 2.11. For all integers n>l, positive numbers a, r, < bj < 1 and an arbitrary 
soe [0,a/2], 

J a<si<-- <.Sn,<a-\-r 

(2.39) 

where > is a constant depending only on a and the bj 's. In particular, if bj = a 
for all j = 1, . . . ,n, then 

p n 

/ n(*^"^^-i)""dsi---ds„<c^_^^(n!)"-V"(i-(i-i/")"). (2.40) 

J a<.si<---<STi<a+r j — i 

Finally, we conclude this section by briefly recalling some aspects of the theory of local 
times. For excellent surveys on local times of random and deterministic vector fields, we 
refer to Geman and Horowitz [20] and Dozzi [17]. 

Let X{t) be a Borel vector field on with values in E''. For any Borel set / C R^, 
the occupation measure of X on / is defined as the following measure on M.'^: 

M.(») = AN{ie/:X(<)G.}. 

If fXj is absolutely continuous with respect to A^, we say that X{t) has local times on 
/ and define its local times, i(», /), as the Radon-Nikodym derivative of fi, with respect 
to A(j, that is, 

In the above, x is the so-called space variable and / is the time variable. Sometimes, we 
write L{x,t) in place of L{x, [0,t]). Note that if X has local times on /, then for every 
Borel set J I, L{x, J) also exists. 

By standard martingale and monotone class arguments, one can deduce that the lo- 
cal times have a measurable modification that satisfies the following occupation density 
formula (see Geman and Horowitz [20], Theorem 6.4): for every Borel set / C and 
every measurable function / : K'' ^ K+ , 



f{X{t))At= / f{x)L{xJ)dx. (2.41) 
Suppose we fix a rectangle / = riiLi[oi,ai + hi]. Then, whenever we can choose a 



version of the local time, still denoted by L{x, Oili [oii + such that it is a continuous 
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function of (a;, ii , . . . , tw) G M'^ x YiiLi [^i ^i] j ^ is said to have a jointly continuous local 
time on /. When a local time is jointly continuous, L{x, •) can be extended to be a finite 
Borel measure supported on the level set 



Xjr\x) = {teI:X{t)=x}; 



(2.42) 



sec Adlcr [1] for details. In other words, local times often act as a natural measure on 
the level sets of X. Hence, they are useful in studying the various fractal properties of 
level sets and inverse images of the vector field X. In this regard, we refer to Berman 
[10], Ehm [19] and Xiao [41]. 

It follows from Gcman and Horowitz [20], ((25.5) and (25.7)) (sec also Pitt [34]) that 
for all x,y E R^, I G A and all integers n>l, 



E[i:(.x,/)"] = (27t) 



7" 



exp {u'' , x) 



and for all even integers n > 2, 

E[(L(.T,/)-L(y,/))"] = (27t)-"'^ 



:cxp^J2{u^,X{t^))j dudt 



(2.43) 



(2.44) 



X EexpJ^i^ {u^,X{t^))j dudt, 



where ?I= (u\.. .,ii"),t= (t\...,r) and each u^ e R'^,P G / C (0,oo)^. In the coordi- 
nate notation, we then write = {u{, . . . , u;^). 



3. Local times: existence and joint continuity 

In this section, we consider the existence and regularity of local times of mfBs. Theorem 
3.1 and Corollary 3.2 provide sufficient conditions for the existence of local times. Then 
in Theorem 3.4 we prove that, under the same condition as in Theorem 3.1, the local 
times of mfBs have a jointly continuous version. 

For I G A, let H — {Hi, . . . , Hn) be the vector defined by 

He = maxHeit) for £ = 1, . . . , TV. (3.1) 
The index H depends on /, but, for simplicity, we have deleted / from the notation. 
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Theorem 3.1. Let I E A and H he the vector defined in (3.1). If d < J^eLi then 
nifBs {B^^^\t)} admits an L^ -integrable local time L{-,I) almost surely. 

Proof. Without loss of generality, we may, and will, assume that S = diam(/) is suffi- 
ciently small that (2.9) holds for all s,t G I. In particular, we assume 6 <l. 

To prove the existence of local time on /, by Geman and Horowitz [20] Theorem 21.9, 
it suffices to prove that 

{E[B^^*\t) - B"'^'\s)fy'^^^dsdt < oo. (3.2) 



I J I 



It follows from (2.9) and a change of variables that 

(E[B^^'\t)-B^^^\s)]Y'^'dsdt<cJj(j^^^^^ ' dsdt 

(3.3) 

(TV \ ^d./'^ 

Since d < y^^-i =-, we can estimate the last integral in a way similar to the proof of 
Xiao and Zhang [45] Theorem 3.6. Namely, by applying Lemma 2.8 repeatedly, we derive 
that the last integral is finite. This proves (3.2) and hence Theorem 3.1. □ 

The following is a consequence of Theorem 3.1 which gives a more natural condition 
for the existence of local times of mfBs. 

Corollary 3.2. Let I e A be fixed. If > d for all t e /, then mfBs (t)} 

admits an L'^ -integrable local time L{-,I) almost surely. 



Proof. Since the functions IIi{t), . . ., H^it) are uniformly continuous on 7, we can di 

1 

dp) 



vide I into subintervals {Ip} such that for each Ip, we have X^fLi -^-^ > d, where 



Heip) = maxtg/p 77£(t). It follows from Theorem 3.1 that on every Ip, mfBs {B^^*\t)} 
has an L^-integrable local time L{-,Ip) almost surely. This implies that {B^'-^^t)} has 
an L^-integrable local time on /, which concludes the proof of Corollary 3.2. □ 



Remark 3.3. The condition in Corollary 3.2 is almost the best possible, in the sense 
that, if d > 'Yl!i=i Hi{t) ^'^^ some t G /, then it can be proven using Geman and Horowitz 
[20] Theorem 21.9 that {S^(*)(<)} has no L^iW^ x fi)-integrable local times on I. In 
the case when d < J2eLi He(t) t G I, but the equality only holds for t in a set 

of Lebesgue measure 0, the existence of local times is rather subtle and requires the 
imposition of further assumptions on (Hi{t), . . . ,HM{t)). Hence, it will not be discussed 
here. 
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We now consider the joint continuity of the local times of {B^^*^{t)}. For conve- 
nience, we first prove that under the same condition as in Theorem 3.1, the local time of 
{5^f(t)(i)} Jias a version that is jointly continuous in both space and time variables. We 
then apply the same argument as in the proof of Corollary 3.2 to show that the same 
conclusion holds, provided X^fci TiJit) ^ t €z I. Our results extend those of Ehm 

[19] for the Brownian sheet and of Ayache, Wu and Xiao [5] for fractional Brownian 
sheets. 



Theorem 3.4. Let {B^^^^t)} = {B"^*\t),te K^} be a multifractional Brownian sheet 

1 

=1 if.' 



with values in R'^. Let L £ A and let H he the vector defined in (3.1). If d < — 



then {B^^*\t)} has a jointly continuous local time on I. 

The main idea for proving Theorem 3.4 is similar to those in Ehm [19], Xiao [41] and 
Ayache, Wu and Xiao [5]. That is, we first apply the Fourier analytic arguments to derive 
estimates for the moments of the local times and then apply a multiparameter version of 
the Kolmogorov continuity theorem (cf. Khoshnevisan [25]). As in Ayache, Wu and Xiao 
[5], the 'one-sided' sectorial local non-deterministic properties of multifractional Liouville 
sheets proved in Section 2 (see Lemma 2.4 and Proposition 2.5) will play important roles 
in deriving moment estimates in Lemmas 3.5 and 3.7 below. However, due to the non- 
stationarity and the lack of two-sided local non-determinism of multifractional Brownian 
sheets, we need to make several modifications in our proofs. 

For convenience, wc further assume 

0<i?i < ••• <i?w < 1. (3.4) 

Lemma 3.5. Assume the conditions of Theorem 3.4 hold. Let r be the unique integer 
in {1, . . . , A^} satisfying 

T.i^'^<ti- (3-5) 

There then exists a positive constant c^.^, depending only on N, d, H and I, such that 
for all X G K.'', all subintervals T = [a,a + (h)] C / with h> small and all integers n>l, 

E[L{x, T)"] < c]^^ {nl f-f^^ h"f^^ , (3.6) 
where (3r = N - t - Hrd + J2''(=i Hr/He- 

Proof. For later use, we will start with an arbitrary closed interval T = Y[^=i['^(:^£ + 
he] C /. It follows from (2.43) and the fact that {Sf (t)}, . . . , {sf ^'^(t)} are indepen- 
dent copies of {B^ ^ (01 ^^lat for all integers n > 1, 



E[L(x,T)"]<(27T)-"^ / III [ exp 



-iVarK:ui<*^)(t^) 



dUk } dt, (3.7) 
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where Uk = (u^, ■ ■ ■ , u^) G K". Fix k = 1, . . . ,d and denote the inner integral in (3.7) by 
jNffe. Then by Lemma 2.6, we have 



/ exp 


N 

i=i 




/ exp 


-l^Var 

^=1 





dUk. 



(3.8) 



Since (3.5) holds, we apply Lemma 2.10 with A = n^^ and q = d to obtain r positive 
numbers pi, . . . ,_Pr > 1 satisfying (2.36) and (2.37). 

Applying the generalized Holder inequality (Hardy [21] page 140), we derive that 



exp 



-|Var(^uiy,(i^) 



\3 = l 



dUk 



cl.li^detCoviY.it'),...^*^)]-' 



/(2pf) 



(3.9) 



e=i 



where the last equality follows from the fact that {Yt{t^), . . . , Yi{t'^)) is a Gaussian vector 
with mean 0. Hence, it follows from (3.7) and (3.9) that 

E[L(x-,r)"]<<_3 f f[[dctCov{Yi{t^),...,Y{f'))]-''^^^P'^dt. (3.10) 

To evaluate the integral in (3.10), we will first integrate [dt\ ■ ■ ■ dt"] for £ — 1, . . . ^t. 
To this end, we use the following fact about multivariate normal distributions: for any 
Gaussian random vector {Zi, . . . , Z„), 



dot Cov(Zi , . . . , Z„) Var(Zi) J| Var(Zj |Zi, . . . , ). 

i=2 



(3.11) 



By the above fact and Lemma 2.4, we can derive that for every £ G {1, . . . , r} and for all 
t\ . . . , i" e T = n^^i [ai, 0,1 + hi] satisfying 



TTf (ri) 



(3.12) 



for some permutation tt^ of {1, . . . , TV}, we have 



det Cov(y,(ti), . . . , y,(r)) > n(^^ 



(3.13) 
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n 

where ij'^*'"'' :=e (recall the decomposition (2.22)). 

We choose £ <\ minja^, 1 < € < N} so that Lemma 2.11 is applicable. It follows from 
(3.12) and (3.13) that 

[detCov(yf(ti), . . . ,y^(i"))]-''/(2w) dt\ • • - dt^ 

<yc"/ TT ^. d4---df^ (3.14) 

In the above, the last inequality follows from (2.40) (recall that Hgd/pi < 1). 
Combining (3.10) and (3.14), we have 

_ T _ JV 

E[L(.T,T)"]<c;\(n!)I^I=i^^''/«[]/i;'^'-('-'/"'^^''/P^^- J| /i,". (3.15) 

We now consider the special case when T = [a, a + (h)], that is, hi = ■ ■ ■ = Jim = h with 
h < S, (3.15) and (2.37) with A = n^^ and q = d together yield 

E[L(x,r)"] < c" („!)ELi^'^/w/i"(^-(i-""')ELi^^''/p^) 

(3.16) 

<c" (n!)^-'^-/i"'^-. 

— .3,9 ^ ' 

This proves (3.6). □ 

Remark 3.6. In the proof of Lemma 3.5, if we apply the generalized Holder inequality 
to the first integral in (3.8) with positive numbers pi, . . . ,pN defined by 

^ Tf 

then the above proof shows that if T C / e ^, then, similarly to (3.15), the following 
inequality holds: 



E[L(x,T)"] < <\^(7i!)^''AAr(r)"(i-'^), (3.17) 



where ly = c?/(X]fci ^) ^ (Oi !)• We will apply this inequality in the proof of Theorem 
3.4 below. 
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Lemma 3.7. Assume the conditions of Theorem 3.4 hold. Let t be the unique integer in 
{1, . . . , N} satisfying (3.5). There then exists a positive and finite constant Cn, depending 
on N, d, H , L and n, such that for all subintervals T = [a, a + (h)] C / with h> small, 
x,y £ R'' with \x — y\ <1, all even integers n>l and all 7 G (0, 1 A 

E[(L(x,T) - L(2/,T))"] < cn\x - ^g^^g^ 

where Ur — Yl\=i ^ ~ 

Proof. Let 7 G (0, 1 A ^) be a constant. Note that by the elementary inequahties 

|e'"- 1| <2i-'^|up foraUweM (3.19) 
and \u + v]"^ < {ul"^ + \vp , we see that for all u^, . . . , u" ,x, y G M'', 



(3.20) 



where the summation ^ is taken over all the sequences {ki, . . . , A:„) € {1, . . . , d}" 
It follows from (2.44), (3.20) and Lemma 2.6 that for every even integer n > 2, 



E[(L(x,T)-L(y,T))"] 

<\x-yrY^ [ dtf n 



KLrcxp 



iVarK](ii-'",B^(*')(t^)) 



du 



, p p n d 



-lVar(^.iBf*^)(t^) 



du 



< \x-y\'' 



< \x-y\'' 



, p p n d 

'E / iiKii'n^^p 



(3.21) 



, P d p n 



iEVar E"i^^(^') 
fci \j=i / 

e=i \j=i / 



du 



du 



dUk 



where 77^ = 1 if fc = kj and rj^, = otherwise. Note that for every j £ {1, . . . , n}, we have 

WetakePi=He (1 < ^ < iV), A = 1/n and g = d in Lemma 2.10 and let (^ = 1,...,t) 
be the constants satisfying (2.36) and (2.37). Observe that since 7 G (0, ^), there exists 
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an io G {1, . . . ,t} such that 



2i/f„7<l. 



(3.22) 



Combining (3.21) with the generahzcd Holder inequahty. we have that 
E[(L(x,T)-L(y,T))"] 



fc=i 



J]|uJpr;iwo cxp 



\3 = ^ 



(3.23) 



xni 

mo 



exp 



iVarK]^.iy,(^^■) 



\ i/pf - 



For any n points t^, . . . , t" G T, let tti, . . . , ttat be TV permutations of {1, 2, . . . , n} such 
that for every 1 <£< N, 



Let 



^^0= / n^r^^oexp 



iVar|^^^.iy,„(<^) 



dC/fc. 



(3.24) 
(3.25) 



By changing the variables of the above integral by means of the transformation 



j = l,...,n; u 



we have that 



where r^o(o) =0. 

Furthermore, by the elementary inequahty that for ^ > 0, \a — b\^ < c^{\a\^ + \b\^), 
where = 2«-i if ^ > 1 and 1 if ^ < 1, we have that 



(3.26) 
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Moreover, the last product is equal to a finite sum of terms, each of the form 
n;=i Hr'^^^o-^ , where e, = 0, 1, or 2 and E;=i ELi vle, - n. 

Let cr^^^^. =E[(Y^^(i'^«oW)) - Y(,Jt'^«oO-i)))2]. gy Proposition 2.5, wc know that M£„ is 
dominated by the sum over all possible choices of (ei, . . . , £„) G {0, 1, 2}" of the following 
terms 

/ nkir''^^'°^^expf-^f]K)V,2„,,')dl4, (3.27) 

where Vk = (w^, . . . e K". By another change of variable, wj, ~ o'iojvj., the integral 
in (3.27) can be represented by 



fi-^r"'"'"'' I ni-ir''^-^^«^^expf-^f](.„i)^')dM/. 



(3.28) 



3 = 1 



where 



-.1 = X„ n IK-r"^"^"'^ exp(^-^^(^i)2j AWu 



Cn,l — 

is a constant depending on n. Thus, we have obtained that 

n 

M,„<c„n<r''^-"°^^'- (3.29) 

The other integrals for I^Iq in (3.23) are easier and can be estimated similarly. 

Combining (3.23) with Lemma 2.4 (which gives lower bounds for o-^^), (3.28), (3.29) 
and the definition of T, we have 

E[(i(z,r)-i(y,r))"] 



n«o i=i 



X 



< c„|x - vrY. \ L n(C^'^ - t;'/'-'^)-'^"'"''/'''"^^^"'"'^dti ■ ■ ■ dti (3.30) 

TT / TT — n r^T^ d4---di^'x n h'}\. 

In the above, Hi = {«£ < tj'*'^'' < • • • < tj*''"'' < + ft.^} for every 1 < £ < r. 
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We now take hi = ■ ■ ■ = = h < 5 in (3.30). Then by Lemma 2.11 and noting that 
(3.22) holds, we have 

E[(L(a;,r)-i(y,r))"] <c„|x-y|"''/i"('^-(^-")SLi^^''/w-^'o7) 

(3.31) 

where the last inequality follows from Lemma 2.10 and the fact that H(„ < 11^. □ 

The proof of Theorem 3.4 is similar to the proofs of Xiao and Zhang [45] Theorem 4.1 
and Ayache, Wu and Xiao [5] Theorem 3.1; we include it here for completeness. 

Proof Theorem 3.4. Let I ^ [a,b] E A be fixed and assume that (3.5) holds. It follows 
from Lemma 3.7 and the multiparameter version of Kolmogorov's continuity theorem 
(cf. Khoshnevisan [25]) that for every TCI, the mfBs {B^^^^} has almost surely a local 
time L{x,T) that is continuous for all x € M''. 

To prove the joint continuity, observe that for all x,y E and s,t E I with |s — 1| < 6, 
where i5 > is the same as in Lemma 2.2, we have 

E[{L{x,[a,s])-L{y,[a,t])r] 

(3 32) 

< 2"-i{E[(i(x, [a, s]) - Lix, [a,f]))"] +E[(L(.T, [a,t]) - L{y, [a,t]))"]}. 

Since the difference i(x, [a, s\) — L{x, \a,t]) can be written as a sum of a finite number 
(which only depends on N) of terms of the form L{x,Tj), where each Tj S ^ is a closed 
subinterval of / with at least one edge length < |s — 1| , we can use Lemma 3.5 and Remark 
3.6 to bound the first term in (3.32). On the other hand, the second term in (3.32) can be 
dealt with using Lemma 3.7, as above. Consequently, for some small 7 E (0, 1), the right- 
hand side of (3.32) is bounded by c"^^ {\x — 7/| + |s — i|)"'', where n > 2 is an arbitrary 
even integer. Therefore, the joint continuity of the local times on / again follows from 
the multiparameter version of Kolmogorov's continuity theorem. This finishes the proof 
of Theorem 3.4. □ 

Similarly to the proof of Corollary 3.2, we derive from Theorem 3.4 the following, more 
general, result. 

Corollary 3.8. Let I E A be fixed. If > d for all t E I, then mfBs {t)} 

has a jointly continuous local time on I almost surely. 

The next corollary is a direct consequence of Corollary 3.8 and the continuity of the 
Hurst functional. We state it to emphasize that the joint continuity of the local time is 
a local property depending on H{t). 

Corollary 3.9. Lett"E[e,A]^ be fixed. IfY.e=iH^>(^' then mfBs {B^'-^^t)} has a 
jointly continuous local time on U{t^,r^)) for some > almost surely, where U(t'^,rQ) 
is the open ball centered at i° with radius tq . 
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Remark 3.10. Let {B^^*\t),t G M.'^} be a real-valued multifractional Brownian sheet 
defined by the harmonizable representation (of. Herbin [22]) 

^o^'Ht) = n |,,|..,).r/2 ^(dA),W e (3.33) 

where W is the complex- valued Gaussian random measure in M.^ with the Lebesgue 
measure Ajv as its control measure. Based on Dobn'c and Ojeda [16] (see also Stoev 

'^Hit) H(t) 

and Taqqu [35]), {Bq ^ {t)} and {B'q ^ {t)} defined by equation (1.4) are equivalent 
up to a multiplicative deterministic function. For {Bq (t)}, we can prove a similar 
result as in Lemma 2.2 (cf. Boufoussi, Dozzi and Guerbaz [13], Wu [38]). Hence, the 
existence of local times of the (A^, (i)-multifractional Brownian sheet {i?^'^*)(t)} defined 
by (1.5), where {B^^*\t)}, . . . , {Bf^*\t)} are independent copies of {B^'*^(t)}, can be 
proven as in Theorem 3.1 and Corollary 3.2. However, we have not been able to prove 
that {Bq ^ (t)} is (even one-sided) sectorial locally non-deterministic (since the Fourier 
analytic technique employed in Wu and Xiao [39] for proving that fractional Brownian 
sheets are sectorial local non-deterministic fails when H{t) varies in t). It is an open 
problem to establish the joint continuity of the local times of {i?^^*'' (t)}. 

4. Holder conditions for L{x, •) 

In this section, we investigate the local and global asymptotic behavior of the local time 
L{x, •) at a; as a random measure on M.^. Results in this section carry information about 
fractal properties of the sample functions of mfBs; see Section 5. 

By applying Lemma 3.5, we can prove the following technical lemmas, which will be 
useful in this section. 

Lemma 4.1. Under the conditions of Lemma 3.5, there exists a positive and finite con- 
stant J , depending only on N, d, H and I , such that for all a ^ I and hypercubes 
T = [a,a + (r)] C / with r < 6, where 6 > is the same as in Lemma 2.2, x € M'' and all 
integers n>l, 

E[L{x + B"<^''\a),Tf] < c"^ ^{nl)^-^-r"''^- . (4.1) 

Proof. The proof is similar to the proof of Ayache, Wu and Xiao [5] Lemma 3.11 and 
we include it here for completeness. For each fixed a € /, we define the Gaussian random 
field Y = {Y{t),te } with values in R'^ by Y{t) = B"^*\t) - B^(")(a). It follows from 
(2.41) that if {_B^'^*)(t)} has a local time L{x,S) on any Borel set S, then Y also has 
a local time L{x,S) on S and, moreover, L{x -\- B^^°-\a), S) = L{x,S). With a little 
modification, the proof of Lemma 3.5 works for the Gaussian field Y. Hence, we derive 
that (4.1) holds. □ 

The following lemma is a consequence of Lemma 4.1 and Chebyshev's inequality. 
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Lemma 4.2. Under the conditions of Lemma 3.5, there exist positive constants c^^, b 
(depending only on N , d, I and H ) such that for all aG I, T = [a, a+ (r)] with r G (0, S), 
where 6 > is the same as in Lemma 2.2, x £ and u> 1 large enough, we have 

F{L{x + B"^''\a),T)>c,y^u^~^^-}<exp{-bu). (4.2) 

Let U{t,r) be the open ball centered at t with radius r, let be the vector defined 
by (3.1) with / = U{t, r) and let t{U) be the positive integer satisfying the corresponding 
condition (3.5). By applying Lemma 3.5 and the Borel-Cantelli lemma, one can easily 
derive the following law of the iterated logarithm for the local time L(x, •): there exists 
a positive constant 3 such that for every x G M'^ and t € (0, 00)^, 

r L{x,U{t,r)) ^ 

hmsup ^ <C43, 4.3) 

where (/^^(r) = r'^-(!;') (loglog(l/r))^"^-('^) with 

r{U) jjU 

P^^U) -N- riU) H%d + ^ 

Because of the continuity of Hi{t) (1 < ^ < N), it can be verified that 

T{U)^T{t) and l3r(u)^Pr[t) asr^O, (4.4) 
where rit) is the unique integer satisfying 

r(t)-l r(t) 



and 

13,^,^=N-T{t)-H,^,){t)d + Y,^f^- (4.6) 

It follows from Fubini's theorem that with probability one, (4.3) holds for almost all 
t € (0, 00)^. We now prove a stronger version of this result, which is useful in determining 
the Hausdorff dimension of the level set. 

Theorem 4.3. Let L E A be a fixed interval and assume that d < J2e=i^- Por any 

fixed X G , let L{x, •) be the local time of {B^^*\t)} at x which is a random measure 
supported on the level set. There then exists a positive and finite constant 4, independent 
of X, such that with probability \, the following holds for L{x, ■)-almost all t G /.■ 



(4.7) 



Local times of multifractional Brownian sheets 889 

where (fi^ir) = r'^^<*) (loglog(l/r))^^^^(') and where (ir(t) defined by (4.6). 

Proof. For every integer fc > 0, we consider the random measure Lk{x,u) on the Borel 
subsets C of / defined by 



L,(.,C)./^(2.trexp(- ^l^""W-^l' )„ 



C JR'^ 



2 



cxp( -^ + i{^,B''W^t)-x)]d^dt. 



(4.8) 



Then, by the occupation density formula (2.41) and the continuity of the function y i— > 
L(?/,C), one can verify that almost surely Lk{x,C) — > L{x,C) as fc — )■ oo for every Borel 
set cel. 

For every integer to > 1, let fm{t) = L{x, U(t, 2""*)). From the proof of Theorem 3.4, 
we can see that almost surely the functions fm(t) are continuous and bounded. Hence, 
we have, almost surely, for all integers m,n> 1, 

/™(t)]"L(x,dt)= hm f[U{t)rLk{xAt). (4.9) 

It follows from (4.9), (4.8) and the proof of Pitt [34], Proposition 3.1 that for every 
positive integer n>l, 

j[f,^{t)fL{xAt) 

^ ^ (n+l)d I- I- I- ( ra+1 



27r 



exp 



i^(.T,uJ) (4.10) 



Eexp|^i^(u^B^(''')(s^))j duds. 



where u = (u\ . . . , 1*"+^) G ]R("+i)'' and s = (t, s\ . . . , s"). Similarly to the proof of (3.6), 
for sufRciently large to, we have that the right-hand side of equation (4.10) is at most 



c 



ds 



h h{t.2-^Y [det Cov(B(f {t),B^^''\s^), Bl^^'"\s"))]<i/^ 

(4.11) 

< c" (,^!)^-^T([/)2-™"/3^(c')^ 

where ^ is a positive finite constant depending only on N, d, H, and /. 

Let 7 > be a constant, the value of which will be determined later. We consider the 
random set 

/™(c.) = {te/:/,„(t)> 7^,(2-™)}. 
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Denote by fi^ the restriction of the random measure L{x, •) to /, that is, ijlu){E) = L{x,EC\ 
I) for all Borel sets £; G M^. We now take n = [log mj , where [yj denotes the integer 
part of y. Then by applying (4.11) and by Stirling's formula, we have 

E/,[/„(i)]"i(x,di) 



[7(^„(2-'")]^ 

C" (77,!)^-'3^(I7)2~™"^t(U) 
<C <C TfL~'^ 



(4.12) 



provided 7 > is chosen large enough, say, 7 > 2 ■ This implies that 

Aic^(/m)^ <oo. 

Therefore, with probability 1 for /Xi^-almost all t G /, we have 

i(a;,C/(t,2-™)) ^ 

By (4.4), we can see that for m sufficiently large, there exists a constant c^^ > such 
that (/3„(2-'") < c^ 5V5,(2-™). Therefore, we have 

i(a;,C/(t,2-™)) 
hmsup ^ ' \. <c, ,. (4.14) 

m — *oo '^f- \^ } 

Finally, for any ?' > small enough, there exists an integer m such that 2"'" < ?- < 2~'"+^ 
and (4.13) is applicable. Since tp^ is increasing near 0, (4.7) follows from (4.13) and a 
monotonicity argument. □ 

Recall that the pointwise Holder exponent of a random field {X{t),t G M^} at a point 
i G is defined by 



axit) =snp< a: lim — =0 



(4.15) 



By Theorem 4.3, and noting that L{-,t) vanishes outside some compact set U (depending 
on w), we have the following corollary. 

Corollary 4.4. For every x G M.'^, the pointwise Holder exponent aL of L{x,t) at t 
satisfies 

aL{t)>l3r(t) a.s. (4.16) 
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5. Hausdorff dimensions of the level sets 

For X e M'*, let {t e (0, oo)^ : (t) ^ x} be the level set of the multifractional 

Brownian sheet In this section, we determine the Hausdorff dimension of T^- 

We remark that the corresponding problem of finding the Hausdorff dimensions of the 
level sets of multifractional Brownian motion has been investigated by Boufoussi et al. 
[13] and the Hausdorff dimensions of the level sets of fractional Brownian sheets were 
studied in Ayache and Xiao [6] . As shown by the following theorem, the fractal structure 
of r^; is much richer than the level sets of multifractional Brownian motion and is locally 
reminiscent of the level sets of a fractional Brownian sheet. 

Theorem 5.1. Let {B^'-^^t)} ^ {B"'•*^t),t eR^} be an {N,d) -multifractional Brow- 
nian sheet with Hurst functionals He{t) (£ = 1, . . . , N ). For any interval L Cz A, let i* G / 
be a point satisfying 



and 

0<Fi(i*)<---<i7jv(t*)<l. 

U Sfci He\f) ^ ^''■'^'^ f'^^ every x £ M^, we have r^. n / = a.s. Lf X^fci Hi^{f) 
then for any x £ M.'^ , with positive probability, 



dimjr,n/)=/3,(t.), (5.1) 



where 



f3r(t') = +N-k- Hk{t*)d, 1 < fc < ivj. 

Remark 5.2. It can be verified that if (4.5) holds for t = t*, then Pr{f) is the same as 
in (4.6) with t replaced by t* . In the special case where H{t) = H is a constant, Pr(t') 
reduces to the form derived in Ayache and Xiao [6] for the Hausdorff dimension of the 
level sets of a fractional Brownian sheet. 

When X^fci Hi^{f) ^ believe that for every a; g R'', Ci I = a.s. However, the 

method of this paper is not sufficient for proving this statement. 

Proof of Theorem 5.1. We will follow the proof of Ayache and Xiao [6] Theorem 5. 
First, we prove 

dim„(r,n/)<min|^^^+iV-fc-iJfe(r)d,l<fc<Ar| a.s. (5.2) 



892 



M. Meerschaert, D. Wu and Y. Xiao 



and r^. n / = a.s. when the right-hand side of (5.2) is negative. 

Without loss of any generaUty, wc may assume that / = [a,a + (h)] and h is small 
so that Lemma 2.2 is applicable. For an integer n > 1, divide the interval / into m,i < 

subrectangles i?„j of side lengths ri"i/^^(**)/i = 1, . . . , iV). Let 0< 
g < 1 be fixed and let k„ j- be the lower-left vertex of Rnj ■ Define p{s, t) = (E[i3^^*'' (t) — 
5ff(s)(s)]2-)i/2^ denote by Np{Rnj,e) the smallest number of p-balls of radius e needed 
to cover Rn.j and denote by D the diameter of Rn.j, that is. 

D := sup{p(s, i):s,t£ Rn,j}- 

By Lemma 2.2, we have 

^ ^\ogNp{Rnj,e)de<cn-\ 
Combining the above inequality with Talagrand [36], Lemma 2.1 gives 

p( max >7i-(i-^)l ^c-™''. (5.3) 

Then for n sufficiently large, we have, for any fixed a; € /, 

P{xe{B"^"Hu),ueRn.j}} 

<p( max \B"^''\s)~B"^''>{t)\<n-'^^-B'>;xe{B"^''\u),ueRn,j}\ 
l^s.teRnj ' J 

+ p( max |S^(^'(s)-B^(*)(t)| >n^(i-e4 (5.4) 
< P{|B^('^-^)(k„j) -x\< + 

M (t) 

In the above, we have applied the fact that Var(_BQ (t)) > c for all t e / to derive the 
last inequality. 

If J2e=i He\f) ^ choose q> such that (1 — g)d > J2i=i He{t') ■ ^'^^ ^'^ 

number of rectangles Rnj such that x£ {B^*^"^ (m), u g Rn,j}- It follows from (5.4) that 

E(^n)<C5,nSr=i(-f^*(**))"'n-(i-e)'^^0 asn^oo. (5.5) 

Since the random variables iV„ are integer-valued, (5.5) and Fatou's lemma imply that 
a.s. Nn = for infinitely many n. Therefore, C] I = a.s. 

We now assume X^fci Hi\f) ^'^^ every fc e {1, 2, . . . , N}, define 
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By Ayache and Xiao [6], Lemma 7, we have 77^. > 0. Define a covering {R'^ j} of F^; n / 
by R'^ J = Rn.j a X £ (m), u g Rn,j} and i?^ ^ = otherwise, i?^ j can be covered 

by n^*='=+i^^'''"* ' cubes of side length ) h. Thus, we can cover the 

level set F^: n / by a sequence of cubes of side length n"-^''^*'^ h. Let g G (0, 1) be small 
and let 



V. 



Hkjt*) 



+ N - k - Hk{t*){l - Q)d. 



Clearly, rj^ [g) > t^^ > 0. In the following, we prove that the Hausdorff dimension of F^; n / 
is bounded above by ?7^(^?) almost surely. To this end, by (5.4), we have 



E 



"ItI at 



(")(M),tiGfl,.,j}} 



(5.6) 



Fatou's lemma implies that the ij^ (g) -dimensional Hausdorff measure of F^. n / is finite a.s. 
and thus dim^ (F^: H I) < rj^, (g) almost surely. Letting g 10 along the rational numbers, 
we obtain dimjj(F2; n /) < rj^ and therefore (5.2). 

To prove the lower bound in (5.1), we assume T{t*) ~ k, that is. 



fc-i 



T.H,{t*)-'^^^^m{t*y 



By Condition A, there exists a positive number such that for all t E [t* — (<;),t* 
k)] n /, we have 



fc-i 



(5.7) 



where Hi (1 < ^ < N) are defined as in (3.1) with in place of /. Note that (5.7) and 
Ayache and Xiao [6], Lemma 3.3, imply that 

y2^+N^k~Hkde{N~k,N^k + l]. 

1=1 



Thus, wc can choose g>Q such that 
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It is sufficient to prove tfiat tfiere is a constant c,^ 3 > such that 

iim,(r,n/,)>7}>c,,3. (5.9) 



Our proof of (5.9) is based on the capacity argument due to Kahane [23]. Similar methods 
have been used by Adler [1], Testard [37] and Xiao [40], to mention just a few. 

Let be the space of all non-negative measures on with finite 7-energy. It is 
known (cf. Adler [1]) that is a complete metric space under the metric 

MU-f I (5.10) 

We define a sequence of random positive measures on the Borel sets C C by 

?H(t)f-f\ _ ^\2- 



M.(C)=X(2.n)^/^cxp(- "l«""'W-l^ )d, 



exp 



(5.11) 

2n 



It follows from Kahane [23] or Testard [37] that if there are positive and finite constants 
C5 4 and C5 5 , independent of g and such that 

E(||a*„||)>c,,„ E(||a*„||2)<c,,3, (5.12) 
E(||/^n||7)<+oo, (5.13) 

where \\fJ.n\\ = /^n(-^?), then there is a subsequence of {/in}, say {/int}, such that fi„k ~* M 
in and /i is strictly positive with probability > ^/{2c^ 5)- It follows from (4.8) that 
^ has its support in n /<; almost surely. Hence, Frostman's theorem yields (5.9). 

It remains to verify (5.12) and (5.13). Let a'^{t) = Var(B^(*)(t)). By Fubini's theorem, 
we have 

E(||/i„il) = y^ ^^e-'<«^->exp(^-|l^)Eexp(i(e,i?^(*)(0»d?dt 
-<«^-)exp(^-l(n-i+a2(i))|e|2^ d^dt 



c 



2n ^ 



cxp -777— r— ^TTTTT 



a^t)J 2(n-i+(T2(t)) 



(5.14) 



Denote by l2d the identity matrix of order 2d, by Cov(i?^(''') (s), S^^^*^ (t)) the covariance 
matrix of {B"^-'\s), B^^^^t)), let = n'^hd + CoviB^^-'^s), B"'^*\t)) and denote by 



27t ^ 
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(^,77)' the transpose of the row vector (^,77). Then 

E(llMnin=/ / / / e"'<«+'''->expf-i(e,77)S(e,77)'')d^d77dsdi 

J I, J I, JW JS.'i \ ^ J 

^J^J^ ^^c^v(^-^{x,x)j:-\x,xY^dsdt (5.15) 

< / / f777> dsdt. 

J I, J I, [det Cov(i?o^(^' (s), 

By applying Lemma 2.4 and the regularity of the H{-), it can be proven that for s, t G I^, 

N ^ _ 

det Cov(i?f > c,,, ^ |,s, - t.l^^^W > c,,, ^ - t,f"'. (5.16) 

Combining (5.7), (5.15), (5.16) and repeatedly applying Lemma 2.8, we obtain 

E(||/x„|l2)<c, A / \ dsdt:^c,^,<oo. (5.17) 

Thus, we have shown that (5.12) holds. 
Similarly to (5.15), we have 

E(ll.J,)^/J^^/^J^^e-.««...„p(-lK,,)E«,,r)d^d, 

<c,8 / / 1^ Kj = -dsdi (5.18) 



1 



where the two inequalities follow from (5.16) and a change of variables. By using Lemma 
2.9 in the same way, we see that E(||/^„||^) < +00 for any 7 defined in (5.8). This proves 
(5.13). 

Finally, by letting g J, 0, the lower bound for the Hausdorff dimension in (5.1) follows 
from (5.9) and we have therefore proven Theorem 5.1. □ 

The proof of Theorem 5.1 suggests that we can consider the Hausdorff dimension of the 
level set in any neighborhood of a point t € (0,oo), provided X^fci TlJit) ^ ^- However, 
in order to obtain an almost sure result, we have to consider at a random level 
X = B^^^\t). The following corollary can be considered as a local Hausdorff dimension 
result for the level sets of mfBs. 
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Corollary 5.3. Let {B"^*\t)} = {B"^*\t),t gM.^} be an {N,d)-multifractwnal Brow- 
ian 

^e=i H^(tO) > d, then there exists vq > such that 

limdim„(rB«(t)(t)n[/(<°,r))=/3,(to) for a.e. tef/(<°,ro) \ = 1, (5.19) 
where /3,-(tO) defined by (4.6) with t'^ in place oft. 



man sheet with Hurst functionals Hg{t) (£ = 1,...,N). If t^ G (0,oo)^ satisfies 



Proof. For any £ (0, oo)^ such that ^^(^0) > d, there exists a positive number rg 

such that for all s G U{t°, tq), we have J^'Li 777(^ > d. By Corollary 3.9, mffls {t)} 
has a jointly continuous local time on U{t'~',rQ). 

By (5.2), we have that for every < r < and x G M'*, 

dim^(T^r]U{t^,r)) < max /3^(^) a.s. (5.20) 
By (5.20) and Fubini's theorem, we see that 

Pidim„(r^n [/(i°,r)) < max /3^f^^, a.e. a; G M'' I 1. (5.21) 

[ seU{t°,r) J 

Since {i?^(*)(t)} has a local time on U{t°,ro), the occupation density formula (2.41) and 
(5.21) together imply that 

p|dimjj(rBff(t)(t) nf7(i°,r)) < max a.e. i G f7(i°, ro) I = 1. (5.22) 

On the other hand, by using an argument similar to the one used in the proof of 
Berman [9] Theorem 2.1 (see also the proof of Xiao [42] Theorem 1.1), we can show that 
for every e > and r G (0, ro) small enough, 



dimjj(rBH(t)(4)n[/(r,r))> max - e a.e. i G t/(r, tq) ^ = 1. (5.23) 

By letting r I and e i along the rational numbers, we see that (5.19) follows from 
(5.22) and (5.23). □ 

Remark 5.4- Corollary 5.3 shows the explicit way in which the fractal properties of the 
multifractional Brownian sheet vary in space. In short, the local Hausdorff dimension is 
derived from the constant parameter formula. Fractional Brownian sheets are essentially 
fractional integrals of Brownian sheets (cf. Benson et al. [8] and Biermc, Meerschaert 
and Schefher [12]). The term "multifractional" indicates that the order of fractional 
integration varies in space. It would be interesting to explore the connection between 
multifractional Brownian sheets and multifractals. For example, are the level sets of the 



Local times of multifractional Brownian sheets 



897 



mfBs multifractals and, if so, how do their structure functions depend on the Hurst index 
function H{t)l 
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